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Abstract 

nterest  in  nonpetroleum  fueled  ground  vehicles  led  Stafford  to 
develop  the  computer  code,  Electric  Vehicle  Simulation  (EVSIM) .  EVSIM 
was  designed  to  predict  the  performance  of  current  electric  vehicles  or 
to  be  used  in  the  design  of  electric/hybrid  vehicles.  Before  EVSIM 
could  be  used,  it  needed  to  be  verified  by  comparing  its  predictions 
to  the  results  of  a  vehicle  test.  It  was  also  desired  to  improve  the 
code's  ability  to  model  several  types  of  Internal  Combustion  Engines 
that  may  be  used  in  hybrid  vehicles.  The  approach  taken  was  to  test 
the  components  of  the  the  hybrid  vehicle  separately  prior  to  testing 
the  whole  vehicle.  This  was  done  to  verify  separate  sections  of  EVSIM 
prior  to  a  systems  run.  The  results  of  these  comparisons  and  the  system 
comparison  are  included  with  recommended  changes  to  EVSIM. 


I.  Introduction 


Interest  in  alternatives  to  conventional  petroleum  fueled  ground  passenger 
vehicles  in  both  the  private  and  public  sectors  led  to  the  development  of  the 
computer  model,  Electric  Vehicle  Simulation  (EVSIM)  by  Stafford  (Ref  21).  EVSIM 
is  capable  of  modeling  hybrid  (electric/internal  combustion)  powered  vehicles. 
While  this  model  was  extensive  it  lacked  modeling  sophistication  and  sufficient 
validation  in  the  areas  of  battery  use  and  the  auxilary  internal  combus  n  (IC) 
engine . 

Background 

EVSIM  can  be  used  to  predict  the  performance  of  an  existing  electric  or 
hybrid  vehicle,  or  aid  in  the  design  of  a  new  vehicle  based  on  the  performance 
desired  (Ref  21).  The  output  of  the  simulation  includes  such  items  as  average 
cost  of  energy  used  per  kilometer  and  the  amount  of  battery  consumed  per  cycle. 
This  information  taken  with  the  estimated  initial  cost  and  maintenance  cost  can 
be  used  to  establish  economic  levels  of  merit  among  various  systems  of  interest. 
The  vehicle  model  includes  the  effects  of  aerodynamic  drag,  rolling  resistance, 
rotational  and  translational  inertia,  and  power  train  efficiencies.  The  battery 
model  used  in  the  simulation  was  based  on  a  fractional  utilization  algorithm 
with  corrections  for  short  term  discharge  effects  for  lead-acid  deep  cycle 
batteries.  The  driving  cycle  selected  was  the  Federal  Urban  Driving  Sequence 
(Ref  5).  Provisions  were  made  in  the  simulation  to  model  an  internal  combustion 
engine/series  motor,  parallel  hybrid  system. 


-* 


Objectives 

The  initial  objective  was  to  select  and  install  an  IC  engine  into  an  ex¬ 
isting  test  vehicle.  The  engine  had  to  be  able  to  provide  the  required  amount 
of  power  and  still  fit  in  the  test  vehicle  with  as  few  modifications  to  the 
engine  or  the  test  vehicle  as  possible. 

The  second  objective  was  to  use  EVSIM  to  predict  the  performance  of  the 
test  vehicle.  These  results  could  then  be  compared  to  the  actual  results 
obtained  from  the  test  vehicle  to  either  validate  the  computer  simulation  or  to 
point  out  the  areas  where  further  refinements  of  EVSIM  were  needed  to  improve 
its  accuracy  as  a  tool  for  the  evaluation  of  hybrid  systems. 

Approach 

To  meet  the  first  objective  EVSIM  was  used  to  determine  the  power  required 
from  the  internal  combustion  engine  during  high  speed  cruise.  The  engine  speed 
required  of  the  IC  engine  was  determined  by  the  reduction  ratio  between  the 
electric  motor  and  the  electric  clutch  on  the  test  vehicle  With  these  con¬ 
straints  an  appropriate  IC  engine  was  selected  for  testing.  The  IC  eLgine  was 
run  on  a  dynamometer  to  determine  its  operating  characteristics.  The  results 
from  these  tests  were  compared  to  the  predictions  of  EVSIM  and  the  manufacturer' 
specifications.  A  vehicle  description  and  discussion  of  the  engine 
selection  is  contained  in  Chapter  two,  the  testing  and  modeling  of  the  internal 
combustion  engine  is  described  in  Chapter  four. 

A  first  step  toward  the  second  goal  of  validating  EVSIM  was  to  test  the 
battery  model  in  the  computer  simulation.  To  do  this,  battery  discharge  tests 
at  both  constant  and  cyclic  discharge  rates  were  run.  The  results  of  these 
tests  were  compared  to  the  results  predicted  by  the  battery  model  in  EVSIM.  The 
results  of  the  tests  and  a  description  of  the  battery  model  are  in  Chapter  three 
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By  testing  the  batteries  and  the  IC  engine  separate  from  the  rest  of  the 
vehicle,  inaccuracies  in  EVSIM  caused  by  the  method  in  which  they  were  modeled 
could  be  separated  from  those  introduced  by  other  parts  of  the  simulation. 

The  IC  engine  and  the  magnetic  clutch  were  installed  in  the  test  vehicle 
and  the  instrumentation  needed  for  the  vehicle  test  was  selected.  With  the  IC 
engine  and  instrumentation  installed,  the  vehicle  was  tested  in  the  all 
electric  and  the  gas/electric  hybrid  mode.  The  results  and  comparisons  of 
these  tests  and  the  performance  predictions  from  EVSIM  are  presented  in 
Chapter  five. 


II.  Vehicle  Design 


In  order  to  select  vehicle  components  the  performance  goals  and  physical 
constraints  on  the  test  vehicle  had  to  be  identified.  These  goals  and  con¬ 
straints  would  serve  to  guide  the  component  selection. 

Performance  Goals 

While  EVSIM  is  flexible  enough  to  model  a  wide  array  of  hybrid  systems,  the 
vehicle  seen  as  the  target  size  for  this  study  was  a  small  four  passenger  sedan 
(Ref  21).  This  type  of  vehicle  is  used  by  the  USAF  for  a  variety  of  on  and 
off  base  missions.  The  performance  goals  listed  in  Table  I  were  established  for 
this  type  of  vehicle  based  on  the  desires  of  USAF  vehicle  managers  (Ref  16)  and 
the  analysis  of  some  off  base  mission  requirements.  The  90  km/hr  top  speed  with 
80  km/hr  cruise  allows  limited  suburban  use.  The  0  to  50  km/hr  acceleration 
requirement  meets  the  demands  of  typical  urban  driving  (Ref  14).  For  a  vehicle 
with  a  mass  of  1200  kg  and  a  frontal  area  of  1.9  m  ,  published  data  (Refs  11,13) 
indicate  that  12  kw  of  motor  power  are  required  to  achieve  the  goal  of  90  km/hr 
and  that  22  kw  of  power  are  needed  to  reach  the  0  to  50  km/hr  acceleration 
requirement.  The  critical  factor  in  electric  vehicle  design  is  the  total  energy 
storage  capacity,  not  the  maximum  power  available  (Ref  8).  The  estimated  power 
requirement  of  10.5  kw  of  power  for  the  80  km/hr  cruise  times  the  two  hour 
cruise  requirement  results  in  an  estimate  of  21  kwh  of  energy  storage  required. 
These  power  and  energy  requirements  became  important  criteria  for  component 


selection. 
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TABLE  I 

Performance  Requirements 

Payload 

4  Passengers  (270  kg) 

Top  Speed 

90  km/hr 

Acceleration 

0  to  50  km/hr 

15  sec 

Range 

Urban 

80  km 

80  km/hr  cruise 

160  km 

Battery  Selection 

While  the  vehicle  had  to  meet  the  performance  requirements  of  the  previous 
section  it  was  also  constrained  to  use  readily  available  components.  This 
constraint  of  availability  ruled  out  a  number  of  battery  types  (such  as 
zinc-chloride  and  sodium-sulfur)  despite  their  potential  high  performance 
(Refs  7,  8,  14,  17)  due  to  their  current  lack  of  development  and  non¬ 
availability. 

In  lead-acid  batteries  the  goals  of  high  energy  density  and  high  power 
density  can  be  achieved  in  a  compromise  simultaneously  (Refs  12,  17).  Standard 
automotive  batteries  have  high  power  densities  (up  to  200  w/kg  during  starting, 
Ref  14)  but  lack  the  ability  to  recover  from  repeated  deep  cycles.  Industrial 
batteries  are  designed  for  use  in  applications  where  size  and  weight  constraints 
are  not  critical  and  the  maximum  power  demand  is  not  much  above  the  average 
demand.  They  can  sustain  up  to  2000  cycles  but  have  power  densities  limited 
to  below  20  w/kg  (Ref  14).  Golfcart  batteries  were  developed  to  provide 
reasonable  power  densities  (up  to  100  w/kg)  and  still  maintain  deep-discharge 
life  expectancies  of  200  to  400  cycles  (Refs  10,  14).  These  batteries  rep¬ 
resent  the  best  trade-off  in  lead-acid  batteries. 
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Both  high  power  density  and  high  energy  density  are,  in  theory,  available 
from  two  nickel  based  batteries.  However  in  recent  tests  nickel-zinc  batteries 


have  been  unable  to  meet  their  expected  performance  and  have  shown  poor  deep- 
cycle  performance  (Ref  8).  The  nickel-iron  battery  has  been  demonstrated  to 
have  excellent  deep  cycle  life  (1000+  cycles)  and  high  power  densities  (130  w/kg) 
with  25  percent  higher  energy  density  than  lead-acid  batteries.  The  problem 
with  these  batteries  is  that  they  have  poor  energy  efficiency  (typically  50 
percent) ,  about  one  third  less  than  a  lead-acid  battery  due  to  excessive  hydrogen 
formation  during  charging.  The  most  reasonable  choice  for  batteries  for  a  near 
term  vehicle  was  the  six  volt  golfcart  battery,  due  to  its  availability  and  lack 
of  disqualifing  characteristics. 

In  determining  the  number  of  batteries  to  use  in  the  vehicle  the  limiting 
constraint  was  the  size  of  the  vehicle  chassis.  The  maximum  number  of  batteries 
that  could  be  installed  in  the  vehicle  and  still  maintain  a  four  passenger 
capability  was  twelve.  Based  on  a  standard  29  kg,  six  volt  battery  this  resulted 
in  a  72  volt,  350  kg  battery  pack.  The  energy  capacity  of  this  system  was  found 
by  multiplying  the  energy  density  of  30  wh/kg  (Ref  14)  by  the  mass  of  350  kg  to 
get  10.4  kwh.  This  result  shows  that  in  order  to  reach  the  two  hour  cruise  goal, 
the  battery  mass  alone  could  not  provide  the  21  kwh  of  energy  needed.  The  IC 
portion  of  the  hybrid  would  have  to  make  up  approximately  half  of  the  required 
energy  used  at  cruise. 
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Hybrid  Selection 

The  inability  of  the  batteries  to  provide  the  total  amount  of  energy 
needed  to  meet  the  goal  of  a  two  hour  cruise  led  to  the  selection  of  a  6  kw 
electric  motor  (Ref  21).  To  provide  the  rest  of  the  power  needed  for  the  80  km/hr 
cruise  an  IC  engine  rated  at  approximately  7  hp  was  needed.  With  the  size 
of  the  engine  determined  there  were  two  possible  arrangements  of  the  IC- 
electric  combination  to  consider. 

The  first  gasoline-electric  combination  considered  was  a  series  hybrid 
system.  In  this  system  the  IC  engine  drives  a  generator  which  supplies 
electricity  to  the  electric  motor  directly  or  to  the  batteries,  depending 
on  the  power  required  by  the  road  load.  The  advantage  of  this  configuration 
is  the  simplicity  of  the  IC  control,  as  the  IC  engine  can  be  allowed  to  operate 
at  its  most  efficient  speed. 

The  disadvantages  of  this  arrangement  are  the  extra  weight  it  requires 
and  low  efficiency.  Assuming  efficiencies  of  85,  98,  70  and  96  percent  for 
the  engine,  controller,  generator,  and  engine /genera tor  coupling  respect¬ 
ively,  it  has  been  calculated  that  only  56  percent  of  the  energy  produced 
by  the  IC  engine  reaches  the  transmission  (Ref  21).  If  the  engine-generator 
pair  are  used  to  recharge  the  lead-acid  batteries  the  efficiency  is  reduced 
by  approximately  another  25  percent  due  to  losses  in  the  recharging  of 
lead-acid  batteries  (Ref  14).  The  reduction  in  efficiency  is  compounded 
by  the  extra  weight  added  in  the  form  of  a  larger  IC  engine  needed  to  make 
up  for  the  losses  in  the  system  and  the  generator  needed  to  produce  the 
electricity. 

The  second  combination  considered  was  a  parallel  hybrid  system.  In 
this  combination  the  IC  engine  is  coupled  to  the  electric  motor  at  the  trans¬ 
mission.  While  the  engine  may  not  always  be  operating  at  its  most  efficient 


speed,  the  losses  in  the  system  are  reduced  to  those  in  the  motor-engine 
coupling. 

The  parallel  arrangement  was  selected  for  the  test  vehicle.  Since  the  electric 
motor  could  apparently  meet  the  acceleration  demands,  the  main  function  of  the 
IC  engine  was  to  meet  the  high  speed  cruise  energy  requirements.  In  addition  to 
providing  the  extra  energy  needed  for  the  high  speed  cruise,  the  proper  use  of 
the  IC  engine  would  reduce  the  motor  current  demand,  which  would  increase  the 
battery  range.  The  coupling  ratio  selected  was  determined  to  allow  high 
speed  cruise  in  a  relatively  low  transmission  gear,  allowing  a  high  electric 
motor  shaft  speed,  as  this  reduces  the  current  demand  in  a  series  motor  (Ref  21). 

To  facilitate  the  addition  of  a  microprocessor  control  in  a  later  study, 
it  was  decided  that  an  electrically  operated  magnetic  clutch  should  be  used 


to  couple  the  engine  and  motor.  A  parallel  shaft  V-belt  drive  was  designed. 
V-belt  designs  have  the  advantage  of  high  efficiency  (96  percent)  and  low 
cost  (Refs  19).  The  clutch  selected  was  rated  up  to  approximately  20  hp 


this  resulted  in: 


RPM 


168x10 . 1x49 . 7 
15.1 


=  5585 


(2) 


With  the  necessary  power  and  speed  range  determined  the  next  choice 
was  what  type  of  engine  to  use.  There  are  four  basic  types  of  engines 
readily  available:  spark  ignition  four  stroke  cycle,  spark  ignition  two  stroke  cycle 
and,  two  and  four  stroke  cycle  compression  ignition  engines.  In  general  the  following 
comparisons  can  be  made  between  spark  ignition  (SI)  and  compression  ignition  (Cl) 
engines : 

(1)  Power  output  per  unit  weight-  The  Cl  engine  generally  weighs  five  to 
twenty  pounds  or  more  per  hp  while  the  SI  engine  in  general  weighs  one  to 
seven  pounds  per  hp. 

(2)  Power  output  per  unit  piston  displacement-  This  factor  can  be  used  to 
roughly  compare  the  size  of  the  two  engines.  Most  high  speed  Cl  engines 
will  deliver  about  .3  hp  per  cubic  inch  of  displacement,  compared  to  .5 
to  .9  hp  per  cubic  inch  for  SI  engines.  For  this  reason  the  SI  engine 
will  tend  to  occupy  less  space  than  a  Cl  engine  of  the  same  hp  rating. 

(3)  Acceleration-  The  Cl  engine  will  produce  the  best  acceleration  due  to 
the  use  of  fuel  injection  in  this  type  of  engine.  The  SI  engine  can  over¬ 
come  some  of  this  disadvantage  through  the  use  of  acceleration  pumps. 

(4)  Reliability-  In  general  the  Cl  engine  is  built  to  stand  rougher  duty 
and  is  rated  well  below  its  maximum  power  output.  This  must  be  weighted 
against  the  SI  engine's  easier  starting,  particularly  in  cold  weather. 

(5)  Fuel  Economy-  The  single  greatest  advantage  of  the  Cl  engine  is  its 
superior  fuel  economy  at  both  full  throttle  (10  to  25  percent  greater)  and 
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at  part  throttle.  The  fuel  used  by  the  Cl  engine  tends  to  be  less  expensive 
which  also  make  the  Cl  engine  less  costly  to  operate. 

(6)  Life  cycle  costs-  The  initial  cost  of  Cl  engines  tends  to  be  greater  than 
the  initial  cost  of  a  comparable  SI  engine,  although  this  may  be  offset  by 
the  longer  life  expectancy  of  the  Cl  engine.  The  maintenance  cost  of  the  Cl 
engine  may  be  slighty  higher  than  the  SI  engine. 

(7)  Operating  speed-  The  SI  engine  tends  to  operate  at  higher  speeds  than  the 
Cl  engine,  3000  to  5000  rpm  for  a  typical  SI  engine  compared  to  1200  to  3000 
rpm  for  a  high  speed  Cl  engine  (Ref  15). 

(8)  Miscellaneous  considerations-  The  fuel  used  by  the  Cl  engine  is  less  volatile 
and  thus  safer  to  use  and  the  Cl  engine  is  better  suited  to  two  cycle  operations. 

The  exhaust  from  the  SI  engine  tends  to  have  a  less  objectionable  odor  but  produces 
more  CO. 

The  distinguishing  feature  of  the  two  stroke  cycle  engine  compared  to  the 
four  stroke  cycle  engine  is  that  there  is  one  power  stroke  for  every  revolution 
of  the  crankshaft  in  the  two  stroke  cycle  engine.  Whereas  there  is  one  power 
stroke  for  every  other  revolution  of  the  crankshaft  in  the  four  stroke  cycle 
engine.  This  would  tend  to  indicate  that  a  two  stroke  cycle  engine  could  produce 
twice  as  much  power  as  the  same  size  four  stroke  cycle  engine.  Unfortunately, 
losses  in  power  and  efficiency  occur  during  the  scavenging  process  when  exhaust 
gases  are  removed  from  the  cylinder.  In  most  SI  two  stroke  cycle  engines  some  of 
the  fresh  fuel  and  air  mixture  is  lost  pushing  out  the  exhaust  gases  reducing  the 
engines  fuel  efficiency.  Spark  ignition  two  stroke  cycle  engines  tend  to  be  limited 
to  small  engines  where  this  fuel  loss  is  not  significant.  In  combustion  ignition 
engines  where  the  fuel  is  injected  into  the  cylinder  just  before  it  is  burned,  fuel 
is  not  lost  as  only  air  is  used  in  the  scavenging  process.  The  two  stroke  cycle 
engine  also  tends  to  be  smaller  and  lighter  than  a  comparable  four  stroke  cycle 


engine.  Two  stroke  cycle  engines  also  tend  to  be  less  expensive  due  to 
their  simpler  design  (ref  18,  9). 


In  addition  to  having  to  provide  the  required  7  hp  at  5  to  6  thousand 
rpm  the  engine  selected  for  the  test  vehicle  had  to  be  able  to  fit  in  the 
limited  space  available  for  it.  The  space  requirement  ruled  out  the  Cl 
engines  that  were  being  considered.  The  choice  left  was  between  several 
two  and  four  stroke  cycle  spark  ignition  engines .  The  two  stroke  cycle 
engines  operated  at  their  rated  power  output  between  5  and  6  thousand  rpm 
while  the  four  stroke  cycle  engines  operated  at  3600  rpm.  The  two  stroke 
cycle  engines  also  had  the  advantage  of  weighting  less  than  15  pounds  while 
the  four  stroke  cycle  engines  weighted  between  50  and  60  pounds.  For  these 
reasons  an  8  hp,  two  stroke  cycle  spark  ignition  engine  was  chosen  for  the 
test  vehicle. 

Vehicle  Description 

The  vehicle  selected  to  test  the  EVSIM  code  was  a  light  four  passenger 
sedan.  It  had  a  frontal  area  of  1.8  meters  square  and  a  fixed  mass  of  325 
kilograms.  The  original  four  speed  transmission  was  used  but  the  gasoline 
engine  and  fuel  tank  were  replaced  by  a  6  kilowatt  electric  motor  and  350 
kilograms  of  batteries.  The  electric  motor  is  controlled  with  a  five  step 
voltage  switch.  The  final  allowable  payload  mass  after  the  conversion  to 
electric  power  was  270  kilograms.  The  vehicle  was  also  instrumented  to 
measure  the  current  demand,  and  RPM  of  the  electric  motor;  the  vehicle 
velocity;  engine  RPM;  and  the  fuel  flow  to  engine.  Provisions  were  made  to 
allow  the  inclusion  of  a  microprocessor  to  control  the  vehicle  and  to 
acquire  data. 


Ill .  Battery  Simulation  and  Evaluation 


The  major  problem  in  predicting  the  range  and  performance  of  electric 
and  hybrid  vehicles  is  the  lack  of  accurate  models  of  the  batteries  (Ref  21). 
The  models  used  need  to  be  improved  and  tested  to  determine  their 
limitations. 


Objective 

The  objective  of  this  effort  was  to  test  the  batteries  used  in  the 
test  vehicle  and  compare  their  performance  to  that  predicted  by  the  battery 
model  used  in  HVSIM.  The  test  results  would  be  used  to  modify  the  battery 
model  if  necessary  and  determine  the  limitations  of  the  battery  model. 


Background 

The  modeling  of  the  batteries  used  in  electric  and  electric  hybrid 
vehicles  is  the  most  difficult  problem  in  obtaining  agreement  between  theory 
and  practice.  Compared  to  other  components  of  the  vehicle  there  is  little 
experimental  data  and  few  analytical  models.  Some  of  the  reasons  for  this 
lack  of  information  include: 

1.  The  lack  of  a  need  for  analytical  models  in  most  past  and  present 

non-electric  vehicle  uses.  ] 

2.  The  electrochemical  nature  of  batteries  is  more  complex  to  model  j 

than  the  electromechanical  nature  of  most  of  the  vehicles'  other  components.  j 

Some  electrochemical  equations  have  been  developed  to  describe  the  internal  j 

phenomena;  the  formidable  task  of  relating  these  to  the  batteries  external  j 

characteristics  still  remain.  I 
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3.  There  is  a  wide  range  of  battery  types  and  even  within  a  single 
type,  there  is  considerable  variation  in  geometry,  plate  thickness,  chemical 


reaction,  and  internal  configuration. 

4.  Battery  performance  depends,  in  general,  upon  the  charge/discharge 
history  of  the  individual  battery  (Ref  22). 

With  these  limitations  in  mind  the  battery  model  used  in  EVSIM  was  designed 
to  represent  lead-acid  batteries  only,  and  to  take  into  account  the  dis¬ 
charge  rates  experienced  by  the  batteries. 

Battery  Simulation 

The  energy  available  from  a  battery  at  constant  load  is  primarily 
a  function  of  the  battery's  discharge  rate.  Figure  3.1  is  the  discharge 
to  capacity  curve  used  to  develop  the  model  used  in  EVSIM.  An  equation 
relating  the  batteries  capacity  to  its  constant  discharge  rate  and  the 
temperature  was  developed  and  presented  by  Stafford  (Ref  21).  This  equation 
was  modified  to  take  into  account  the  effects  of  varying  discharge  rates. 
Combining  this  equation  with  one  for  the  internal  resistance  of  the 
batteries  allows  the  terminal  voltage  of  the  battery  to  be  estimated  (Ref  21). 

Model  Verification 

Five  goals  for  the  test  were  established  to  meet  the  objective  of  ver¬ 
ifying  the  battery  simulation  and  making  any  corrections  needed.  The  first 
goal  was  to  verify  that  the  simulation  could  accurately  model  a  battery 
being  discharged  at  a  constant  rate.  The  second  goal  was  to  determine  the 
effect  that  a  temperature  change  has  on  the  simulations  accuracy.  The 
third  goal  was  to  test  the  battery  at  varying  discharge  rates  and  compare 
the  measured  results  to  those  predicted  by  the  computer.  The  fourth  goal 


was  to  measure  the  effect  of  the  battery's  self  regeneration  during  a 
period  of  rest  on  its  useful  energy.  The  fifth  goal  was  to  determine  the 
effect  of  overcharging  a  battery  on  its  range. 

A  schematic  of  the  test  setup  used  is  shown  in  figure  3.2.  For 
the  first  series  of  tests  to  verify  the  sir-ulation  at  a  constant  discharge 
rate  the  battery  was  discharged  through  constant  resistor  Rl.  The  third 
goal  was  achieved  by  using  a  series  of  relays  to  change  the  value  of  Rl 
thereby  varying  the  discharge  rate.  These  tests  were  performed  at  21°C  and 
at  6°C  to  determine  the  effect  of  temperature  on  the  simulation  accuracies.  T 
battery  was  overcharged  by  .3  volts  to  measure  its  effect  on  the  useful  energy 

Results 

From  the  results  of  the  constant  discharge  tests,  it  could  be  seen 
that  the  equations  for  constant  discharge  needed  to  be  modified.  Using 
Stafford's  original  data,  a  new  set  of  equations  were  developed.  When 
the  new  equations  were  used,  the  difference  between  the  measured  data 
and  the  predicted  data  fell  from  over  1  volt  to  less  than  .4  volts.  Using 
these  new  equations  the  rest  of  the  experimental  data  was  compared  to  the 
battery  simulation  predictions.  Table  II  lists  the  results  of  these 
comparisons.  Figures  3.3  and  3.4  show  results  of  typical  constant  rate 
discharge  and  a  varying  rate  discharge  respectively.  Figures  3.5,  3.6,  and 
3.7  show  the  effects  of  discharge  rate,  temperature,  and  overcharging  of 
the  battery  on  the  prediction  accuracy.  Figures  3.8  and  3.9  show  the 
results  if  the  simulation  is  based  on  discharge  equations  generated 
from  the  experimental  data  of  test  one.  A  complete  description  of  the  tests 
and  their  results  are  in  Appendix  A. 
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Figure  3.2  Battery  Test  Schematic 

Conclusions 

1.  The  battery  model  used  in  HVSIM  can  predict  constant  discharge 
battery  performance  within  10%  of  measured  performance.  The  model 
tends  to  over  predict  the  available  energy. 

2.  The  battery  model  predicts  varying  discharge  rate  performance 
within  10%. 

3.  Using  the  battery  model  based  on  the  measured  battery  performance, 
the  difference  between  measured  and  predicted  performance  was  within  5% 
for  constant  or  varying  rate  discharge. 

4.  The  air  temperature  did  not  significantly  affect  the  accuracy 
of  the  battery  performance  prediction. 

5.  The  battery  model  was  able  to  predict  the  effects  of  battery 
regeneration;  however,  regeneration  will  not  significantly  increase  the 
range  of  hybrid  vehicles. 

6.  The  general  voltage  expression  for  a  nominal  six-volt  battery 
given  in  EVSIM  (Ref  21)  should  be  replace  with: 

Vb  =  6.09  -  (.003)Rb*Ib  -  ( . 36)DF2 
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BATTERY  TEST  POUR 


Figure  3.4 


Varying  Discharge  Rate 


Figure  3.6 
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BATTERY  TEST  POUR 


Figure  3.9 


Varying  Discharge  Rate 


IV.  Hybrid  Simulation 

While  a  simple  model  for  an  internal  combustion  engine  was  included 
in  the  original  version  of  EVSIM,  a  more  accurate  and  versatile  model  was 
needed  to  predict  the  performance  of  a  hybrid  vehicle  in  HVSIM. 

Objectives 

The  first  objective  was  to  develop  a  model  that  could  predict  both 
the  available  power  and  the  fuel  consumption  of  the  internal  combustion 
engine.  The  model  had  to  be  versatile  enough  to  cover  both  spark  ignition 
and  compression  ignition  engines. 

The  second  objective  was  to  test  the  internal  combustion  engine 
selected  for  the  test  vehicle  and  compare  the  results  of  the  tests  to  the 
predictions  of  the  model. 

Internal  Combustion  Engine  Model 

The  key  to  this  model  is  the  assumption  that  the  full  throttle  torque 
curve  can  be  represented  by  a  second  degree  equation  of  the  form: 

T  =  a  +  bN  +  cN2  (3) 

where  T  is  torque  and  N  is  engine  rpm.  This  assumption  is  reasonable  in 

the  normal  operating  region  of  most  engines  (Ref  3).  This  method  can  be  used 

for  both  spark  ignition  and  combustion  ignition  engines,  although  their 

torque  curves  are  different  (fig. 4. 1,4.2) .  Spark  ignition  engines  a,  b,  and  c 

are  found  knowing  dT/dN  =  0  at  N„,  T.  =  T,.  at  N, ,  and  T„  =  T 

6  2 ’  1  bhpmax  1 ’  2  max 

at  Nr 
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For  combustion  ignition  engines,  equation  3  can  be  solved  knowing 


dT/dN  =  0  at  N„,  T,  =  T  at  L,  and  T,  =  T  at  N,.  Using  these 

2*  2  max  2’  1  maxrpm  1  ° 

boundary  conditions  to  solve  equation  3  yields 


-  T2)/(N2  -  2N1N2  +  Nj  ) 

(4) 

b  =  -2CNJ 

(5) 

2 

3  =  T2  +  cN2 

(6) 

Having  solved  equation  3,  the  brake  horsepower  can  be  obtained  from: 


T  x  N  (7) 

5252 


Figure  4.3  shows  the  torque  calculated  from  equation  3  plotted  with 
the  torque  curve  supplied  by  the  IC  engine  manufacturer.  Figure  4.4  shows 
the  BHP  predicted  by  the  engine  model  and  the  manufacturer's  horsepower 
curve. 


IC  Engine  Test 

The  engine  used  in  the  test  vehicle  was  tested  in  accordance  with 
SAE  engine  test  code  SAE  J816b  (Ref  20).  The  engine's  torque  was  measured 
with  a  water  brake  dynamometer  and  its  rpm  with  a  stroboscope.  The 
measured  horsepower  is  plotted  with  the  manufacturer's  estimated  horsepower 
in  figure  4.5.  The  measured  torque  was  used  to  make  a  new  estimate  of 
the  engine's  horsepower;  this  is  plotted  in  figure  4.6. 
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Figure  4.5 
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Figure  4.6 


Conclusions 


1.  The  results  in  figure  4.3  and  4.4  shows  the  importance  of  testing 
the  engine  in  the  configuration  it  will  be  used  in.  The  manufacturer's 
horsepower  curve  was  based  on  an  er  ,ine  with  no  muffling.  The  engine  was 
tested  with  the  exhaust  pipe  and  muffler  used  in  the  test  vehicle. 

2.  The  IC  engine  model  is  accurate  to  within  several  percent  over 
the  range  of  interest  for  full  throttle  operations. 


V.  Hybrid  Vehicle  Test  and  Simulation 


In  this  chapter  the  results  of  complete  vehicle  simulation  are 
compared  to  the  measured  performance  of  the  test  vehicle.  The  simulation 
results  and  the  test  vehicle  performance  will  be  compared  for  hybrid  and 
all  electric  configurations. 

Vehicle  Simulation 

The  program  HVSIM  was  run  using  both  the  general  battery  model 
and  the  model  developed  specifically  for  the  batteries  used  in  the 
test  vehicle  both  of  these  models  use  the  change  recommended  in  chapter  three. 
By  comparing  the  results  of  these  simulations  the  effect  of  the  battery 
model's  accuracy  on  the  overall  simulation  accuracy  could  be  determined. 

The  driving  cycle  used  was  the  test  vehicle's  speed  schedule. 

Vehicle  Test 

The  vehicle  was  driven  around  a  5.2  mile  course.  The  test  course 
contained  approximately  2.5  miles  of  hills  up  to  a  10%  grade.  The  speed 
and  electric  energy  used  were  recorded.  The  vehicle  was  driven  around 
the  course  first  with  the  IC  engine  operating  and  then  with  only  the 
electric  engine.  Table  III  shows  the  speed  schedule  of  the  test  vehicle 
for  both  the  hybrid  and  electric  tests. 

Vehicle  Test  and  Simulation  Results 

Table  IV  shows  the  results  of  the  all  electric  test  and  simulation. 

The  first  simulation  used  the  general  battery  model  developed  in  Chapter  3. 

The  second  simulation  uses  the  battery  model  developed  specifically  for  the 


Table  III 


Vehicle  Speed  Schedule 

Time 

Distance 

Speed 

(minutes) (sec) 

(miles) 

(mph) 

51 

.5 

35 

2 

12 

1.4 

40 

3 

02 

2.0 

43 

4 

03 

2.6 

35 

4 

54 

3.1 

35 

6 

15 

4.0 

40 

7 

05 

4.6 

43 

8 

07 

5.2 

35 

Table  IV 

Electric  Test  and  Simulation  Results 

Test  Vehicle 

HVSIM  Experimental  Battery  Model 

Time  (min) : 

8.1 

8.1 

8.1 

Distance  (km) : 

8.4 

8.2 

8.2 

Energy  used  (kw-hr) 

Electric  Energy: 

1.06 

.8 

.8 

Table  V 

Hybrid  Test  and  Simulation  Results 

Test 

Vehicle 

HVISM 

Time  (min) : 

8.1 

8.1 

Distance  (km) : 

8.4 

8.2 

Energy  used  (kw-hr) 

Electric  Energy: 

.65 

.46 

IC  Engine: 

.41 

.33 

Total  Energy: 

1.06 

.79 

batteries  used  in  the  test  vehicle.  Table  V  shows  the  results  of  the 
hybrid  test  and  simulation. 

Conclusions 

1.  The  difference  in  the  accuracy  of  the  battery  models  used 
had  very  little  effect  on  the  overall  accuracy  of  the  vehicle  simu¬ 
lation. 

2.  The  hybrid  simulation  can  accurately  predict  the  relative 
amount  of  energy  supplied  by  the  electric  motor  and  the  IC  engine. 

3.  The  simulation  doesn't  take  hills  into  account  in  its  prediction 
of  vehicle  performance.  This  accounts  for  most  of  the  underprediction 

of  the  energy  used  by  the  test  vehicle. 

Recommendations 

1.  A  more  accurate  evaluation  of  HVSIM  would  require  further 
testing  of  the  vehicle  on  a  flat  test  track  and  added  instrumentation 
to  measure  the  fuel  flow  to  the  engine. 

2.  The  energy  supplied  by  the  IC  engine  could  be  increased  if  the 
pulley  size  was  adjusted  to  allow  it  to  operate  at  more  efficent  speeds. 
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APPENDIX  A 


Battery  Testing 

The  batteries  were  tested  under  the  following  eight  different 
combinations  of  temperature  and  discharge  rate: 

O 

Test  1:  21  C,  158  amps  constant  discharge  for  15  minutes. 

o 

Test  2:  21  C,  93  amps  constant  discharge  for  30  minutes. 

o 

Test  3:  26  C,  44  amps  constant  discharge  for  one  hour  with  a  six  hour 

rest  and  a  second  discharge  for  17  minutes, 
o 

Test  4:  26  C,  varying  discharge  between  46  and  75  amps  for  24  minutes. 

o 

Test  5:  6  C,  40  amps  constant  discharge  for  40  minutes. 

o 

Test  6:  5  C,  90  amps  constant  discharge  for  20  minutes. 

o 

_  Test  7:  7  C,  varying  discharge  between  103  and  132  amps  for  13  minutes. 

(2 

Test  8:  5  C,  67  constant  discharge  for  49  minutes. 

Figures  A1  through  A8  are  the  results  of  these  tests. 

Figure  A. 9  shows  the  test  setup  used  in  the  discharge  experiments. 
The  batteries  were  discharged  through  a  4  inch  by  4  inch  bar  of  graphite 
3  feet  long.  Terminal  posts  were  placed  along  the  length  of  the  bar.  By 
changing  the  terminals  the  battery  was  connected  to  the  discharge  rate  of 
the  battery  was  adjusted.  The  terminal  voltage  and  the  current  flowing 
through  the  rod  were  recorded  on  a  strip  chart  recorder.  The  data  was 
transferred  to  a  VAX  11/780  with  an  HP  digitizer  for  processing. 
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APPENDIX  B 


HVSIM 

This  appendix  is  a  listing  of  HVSIM  with  the  major  changes  frown 
EVSIM  marked  with  an  *  and  a  sample  of  its  output. 

C  PROGRAM  EVSIM( INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT) 

DIMENSION  VX( 1400) , GR( 10 ) , CSTP ( 5 ) , MSTP( 10 ) , ETAM( 10 , 30 ) , ETAT ( 10, 20) 
C, TSTP ( 10 ) , RBAT ( 5 ) , RBAT1 ( 5 ) ,V(1400) 

REAL  MM, M, MS, KMOT, MSTP, Ml , MCF, MSH 

OPEN  (UNIT=1 , FILE*® 1 EVS . DAT  1 , STATUS= 'OLD ', READONLY) 

OPEN  (UNIT=5 , FILE*® ' CARTEST . DAT ' , STATUS= 1  OLD ' , READONLY) 

900  FORMAT ( 10F5. 1) 

901  FORMAT ( 1 8F4. 1) 

WRITE (6, 902) 

902  FORMAT ( 1H 1, 31X6 7H '  *  'ELECTRIC  VEHICLE  SIMU 

C  L  A  T  I  O  N  '  '  '  ) 

ENTER  VEHICLE  CONSTANTS.  DRAG  COEFF( CD) , FRONTAL  AREA(FA),  PAYLOAD 
MASS (PM) ,  FIXED  MASS (FM) . 

READ ( 1 , * )CD, FA, PM , FM 
WRITE (6, 903) 

903  FORMAT (///20X22HVEHICLE  CONSTANTS  ARE:) 

WRITE (6 , 904 ) CD, FA, PM, FM 

904  FORMAT (3 5X1 7HDRAG  COEFFICIENT*®, F5. 2, /35X13HFRONTAL  AREA®*, F5 . 2, 14H 
CSQUARE  METERS,/ 3 5X13H PAYLOAD  MASS=, F4 . 0, 10H  KILOGRAMS, /3 5X11HFIXED 
C  MASS=,F4.0) 

ENTER  ENVIRONMENTAL  CONSTANTS.  AMBIENT  TEMPERATURE (TMP ) ,  ELECTRIC 
POWER  COST(ECOST) ,  PETROLEUM  FUEL  COST(FCOST). 

READ ( 1 , * )TMP , ECOST, FCOST 
WRITE (6, 905) 

905  FORMAT (/2 0X2 8H ENVIRONMENTAL  CONSTANTS  ARE:) 

WRITE (6, 906) TMP, ECOST, FCOST 

906  FORMAT ( 3 5X 2 0H AMBIENT  TEMPERATURE=, F4 . 0, 8H  CELSIUS , /3 5X23HELECTRIC 
CENERGY  COST=  $,F4.3,9H  PER  KWHR, / 3 5X22HPETROLEUM  FUEL  COST=  $,F4.2 
C, 10H  PER  LITER) 

ENTER  VEHICLE  POWER  REDUCTION  VARIABLES.  ROLLING  RADIUS (RR),  TIRE 
ROLLING  COEFFICIENT (CR) ,  AXLE  RATIO (AR) ,  NUMBER  OF  FORWARD  SPEEDS 
)0  IF  CVT  (IS),  MAX  VEHICLE  DESIGN  SPEED(VMAX) . 

READ(1 , * )RR, CR, AR, IS, VMAX 
WRITE (6, 907 ) 

907  FORMAT (/20X 3 0HPOWER  REDUCTION  VARIABLES  ARE:) 

WRITE (6 , 908 )RR, CR, AR, IS, VMAX 

908  FORMAT ( 3 5X20HTIRE  ROLLING  RADIUS  =  , F6 . 3, 7H  METERS, / 3  5X2 0H ROLLING  CO 
CEFFICIENT=, F7 . 4, /3 5X11HAXLE  RATIO=, F5 . 2 , /3 5X2 1HNUMBER  FORWARD  GEAR 
CS=,  I2,/3  5X17HMAXIMUM  VELOCITY*®,  F4. 0,  6H  KM/ HR) 

IF( IS . EQ. 0)GO  TO  20 
10  CONTINUE 


C  ENTER  TRANSMISSION  RATIOS)lF  MULTISPEED  (GR(IS))  ON  FIRST  CARD.  ON 

C  NEXT  CARD  ENTER  DESIRED  MOTOR  UPSHIFT  AND  DOWNSHIFT  RPM(UPSHFT, 

C  DNSHFT ) . 

C 

READ ( 1 , * ) (GR ( I ) , 1=1 , IS ) 

READ ( 1 , * )UPSHFT, DNSHFT 
WRITE (6 ,909) 

909  FORMAT (/20X39HF I XED  SPEED  TRANSMISSION  VARIABLES  ARE:) 

WRITE (6, 910 )UPSHFT, DNSHFT, GR(1 ) 

910  FORMAT ( 3 5X20HMOTOR  UPSHIFT  SPEED=, F5 . 0, 4H  RPM, /3  5X22HMOTOR  DOWNSHI 
CFT  SPEED=,F5.0,4H  RPM , /3 5X1 2HGEAR  RATIOS=, F6 . 3 ) 

IF  (IS.EQ.l)GO  TO  30 

11  DO  12  1=2, IS 

12  WRITE (6,911 )GR ( I ) 

911  FORMAT (47X,F6.3) 

GO  TO  30 

20  CONTINUE 
C 

C  ENTER  CVT  VARIABLES.  RATED  INPUT  POWER (PRATT)  AND  NUMBER  OF  SPEED 

C  RATIO  STEPS  FOR  INPUT  DATA  ON  FIRST  CARD,  EFFICIENCY  DATA(ETAT) 

C  ON  NEXT  CARD  SERIES:  FIRST  CARD-  SPEED  RATIO  STEPS(TSTP),  FOL- 

C  LOWING  CARDS-  ETAT ) . GT . 0  AND  .LT.l  IN  10F5.1  FOR  EACH  5(  OF  PRATT 

C  FROM  5(  TO100( .  2  CARDS  REQD  FOR  EACH  TSTP  . 

C 

READ(1 , * ) PRATT, ITSTP 

READ ( 1 , * ) (TSTP ( I ) , 1=1 , ITSTP ) 

DO  21  1=1, ITSTP 
DO  21  J=1 , 2 

21  READ(1 , 900 ) (ETAT ( I , ( J-l ) *10+K ) , K=1 , 10 ) 

WRITE (6, 91 2) 

912  FORMAT (/20X49HCONTINUOUSLY  VARIABLE  TRANSMISSION  VARIABLES  ARE:) 
WRITE (6, 913) PRATT 

913  FORMAT ( 3 5X13H POWER  RATING*, F5 . 1 , 3H  KW) 

30  CONTINUE 
C 

C  ENTER  MOTOR  VARIABLES.  MOTOR  TYPE(MTYP) )0=SHUNT, 1=SERIES  ,  K  FAC- 

C  TOR(KMOT),  ARMATURE  RESISTANCE(RARM) ,  FIELD  RESISTANCE(RFLD) , 

C  RATED  POWER ( PRATM ) ,  MAX  POW  MULT(IPMAX),  MAX  CUR(AMAX),  BASE  SPEED 

C  (BSPD),  AND  NUMBER  OF  SPEED  STEPS  FOR  INPUT  DATA  ON  FIRST  CARD. 

C  ENTER  EFFICIENCY (ETAM)  ON  NEXT  CARD  SERIES:  FIRST  CARD-  MOTOR 

C  SPEED  STEP  )RPM  (MSTP),  FOLLOWING  CARDS-  ETAM) . GT . 0  AND 

C  .LT.l  IN  10F5.1  FOR  EACH  10(  OF  PRATM  UP  TO  300()3  CARDS  REQD 

C  FOR  EACH  MSTP  . 

C 

READ ( 1 , * ) MTYP , KMOT, RARM , RFLD , PRATM , IPMAX, AMAX, BSPD, IMSTP 
READ ( 1 , * ) ( MSTP ( I ) , 1=1, IMSTP) 

DO  31  1=1, IMSTP 
DO  31  J=1 , 3 

31  READ ( 1 , 900 ) (ETAM( I , ( J-l ) *1 0+K ) , K=1 , 10 ) 

WRITE (6,914) 

914  FORMAT (/20X 2 0HMOTOR  VARIABLES  ARE:) 

IF ( MTYP . NE . 0 ) GO  TO  33 

32  CONTINUE 
WRITE (6 , 915) 

915  FORMAT (3 5X1 7HMOTOR  TYPE=  SHUNT) 


GO  TO  34 

33  WRITE (6 , 916 ) 

916  FORMAT ( 3 5X18HMOTOR  TYPE=  SERIES) 

34  WRITE (6 , 91 7 )PRATM , BSPD, AMAX 

917  FORMAT (3 5X1 2H RATED  POWER=, F5 . 1 , 3H  KW , / 3 5X1 1HBASE  SPEED=, F5 . 0, 4H  RP 
CM, / 35X21HMAX  ARMATURE  CURRENT= , F 5 . 0 , 5H  AMPS) 

40  CONTINUE 

ENTER  CONTROLLER  VARIABLES.  ENTER  TYPE ( CTYP )  ON  FIRST  CARD:  1  IF 
FINITE  STEP  SWITCHING,  2  IF  CHOPPER. 

READ (1,*) CTYP 

IF ( CTYP . NE . 1 . ) GO  TO  50 

41  CONTINUE 

ENTER  FINITE  STEP  CONTROLLER  VARIABLES,  FIRST  CARD-  NUMBER  OF  VOLT 
AGE  STEPS ( ICSTP ) .  NEXT  CARD-  FRACTION  OF  FULL  VOLTAGE  FOR  EACH 
STEP(CSTP) ,  AND  EFFECTIVE  SOURCE  RESISTANCE ( RBAT )  FOR  EACH  STEP. 

READ (1,*) ICSTP 

READ(1,*) (CSTP(I ) ,1=1, ICSTP) , (RBAT (I ) , 1=1, ICSTP) 

WRITE (6,920) ICSTP, CSTP (1 ) , RBAT ( 1 ) 

DO  42  1=2, ICSTP 

42  WRITE ( 6 , 92 1 ) C  STP ( I ) , RBAT ( I ) 

920  FORMAT (/2 0X1 2, 3 1H  STEP  CONTROLLER  VARIABLES  ARE :, /3 5X1 7H VOLT AGE  FR 
CACTION=,F5. 3, 3X17HSOURCE  IMPEDANCE=, F6. 4, 5H  OHMS) 

921  FORMAT ( 5  2XF  5.3, 20XF  5.3) 

GO  TO  60 

50  CONTINUE 

ENTER  CHOPPER  CONTROLLER  SOURCE  IMPEDANCE (RBAT) .  MODEL  ASSUMES  1.5 
VOLTS  JUNCTION  LOSS  AND  3(  RESISTIVE  LOSS. 

READ ( 1 , * ) RBAT ( 1 ) 

WRITE (6,929) RBAT ( 1 ) 

929  FORMAT (/2 0X3 3HCHOPPER  CONTROLLER  VARIABLES  ARE : , / 3 5X1 7HSOURCE  IMPE 
CDANCE=,F5. 3, 5H  OHMS, /3 5X, 32HASSUMED  JUNCTION  DROP=  1.5  VOLTS, /35X2 
C8HASSUMED  RESISTIVE  LOSSES=  3 ( ) 

CSTP(1 )=1 . 

ICSTP=1 

60  DO  61  1=1, ICSTP 

61  RBAT1 ( I )=RBAT ( I ) 

ENTER  BATTERY  VARIABLES.  FIRST  CARD-  NOMINAL  TOTAL  VOLTAGE (VNOM) 
AND  TOTAL  MASS(BM) . 

READ ( 1, *) VNOM, BM 
BM1=BM 

WRITE ( 6, 92 6) VNOM, BM 

926  FORMAT (/2 0X2 2HBATTERY  VARIABLES  ARE : , /3 5X, 16HNOMINAL  VOLTAGE=, F4 . 0 
C, /35X13HBATTERY  MASS=,F5.0) 

ENTER  HYBRID  VARIABLES.  FIRST  CARD-  1  IF  HYBRID, 0  IF  NOT.  NEXT  CAR 
) IF  IHYB=1  -  HYBRID  SPEED  RATIO(HR) ,  COUPLING  EFFICIENCY(ETAH) , 
ENGAGEMENT  SPEED(VH),  NORMAL  OPERATING  TRANSMISSION  GEAR(ISH) , 


C  MAX  TORQUE (TORI ) ,  TORQUE  AT  MAX  RPM(TOR2),  RPM  AT  MAX  TORQUE (EN1) 

C  MAX  RPM( EN2 ) ,  FUEL  FLOW  AT  MAX  TORQUE (FF1),  AND  FUEL  FLOW  AT  MAX 

C  RPM( FF2 ) . 

C 

READ ( 1 ,  *  )  I H YB 
HM=0 . 

VH=100. 

IF ( IHYB . EQ. 0)GO  TO  71 

READ ( 1 ,  * )HR, ETAH, VH, ISH,HM 

READ ( 1 , * )TORl , TOR2 , EN1 , EN2 , FF1 , FF2 

WRITE ( 6 , 92  7 ) HR, ETAH, VH , ISH, HM , TORI , TOR2 , EN1 , EN2 , FF1 , FF2 
927  FORMAT ( /20X, 34H7KW  PARALLEL  HYBRID  VARIABLES  ARE : , /3 5X25HHYBRID/MO 
CTOR  SPEED  RATIO=,F6. 3, /3 5X20HCOUPLING  EFFICIENCY*1,  F5 . 1 ,  1H%,  ,/3 5X17H 
CENGAGEMENT  SPEED**,  F5 . 1 ,  5HKM/ HR,  / 3  5X22HTRANSMISSION  POSITION=,  12 , 

C'  GEAR' ,/35X, ’HYBRID  MASS= '  ,  F5 . 1 ,  3H  KG,/35X,  ’MAX  TORQUE**  '  ,  F5 . 2, /3 5 
CX, 'TORQUE  AT  MAX  RPM= '  ,  F5 . 2 ,  / 3  5X,  '  RPM  FOR  MAX  TORQUE**  '  ,  F7 . 1 , /3 5X,  ' 
CMAX  RPM**' ,F7.  1,/35X, 'FUEL  FLOW  AT  MAX  TORQUE** '  , /3  5X,  '  FUEL  FLOW  AT 
CMAX  RPM= ’ ) 

COEF3=(T0R2-TOR1)/(EN1**2-2.*EN1*EN2+EN2**2) 

COEF2=-2.*COEF3*ENl 
C0EF1=T0R1+C0EF3*EN1**2 
70  IS=*ISH 
'  71  CONTINUE 

i 

!  PROGRAM  IS  NOW  READING  INPUTTED  DRIVING  CYCLE  SPEED  SCHEDULE  DATA. 

:  FIRST  CARD-  TIME  INCREMENT (T ) ,  NUMBER  OF  DATA  POINTS (NDATA) .  FOL- 

!  LOWING  CARDS-  SPEED  SCHEDULE(VX)  IN  18F4.1.  LAST  CARD-  CONVERSION 

!  FACTOR ) 0  IF  VX  IN  MPH,  1  IF  VX  IN  KPH  (CONV) . 

t 

READ(1,*)T, NDATA 
NNDATA*=NDATA/ 1 8+1 
DO  80  J=1 . NNDATA 

80  READ (5 , 901 ) (VX{ 18* ( J-l )+K ) , K=2 ,19) 

READ (1,*) CONV 

PRINT*,  'CONV', CONV 
IF ( CONV . NE . 0 . ) GO  TO  82 

81  DO  8101  J-l, NDATA 
8101  VX(J+1)**VX(J+1)  *1.602 

82  NNNDATA**NDATA/60+l 

ENTER  PROGRAM  OPTION.  FIRST  CARD-  0  IF  SYSTEM  EVALUATION,  1  IF 
SYSTEM  DESIGN.  NEXT  CARD) IF  IPROG-1  -  MAXIMUM  ALLOWABLE  BATTERY 
DISCHARGE ( DSMAX ) ,  AND  MAXIMUM  ALLOWABLE  SPEED  SCHEDULE  DEVIATION 
(DVMAX) . 

READ ( 1 , * ) I PROG 

IF( IPROG. NE . 0 )GO  TO  8202 

8201  WRITE (6, 930) 

930  FORMAT (/2 0X3 3HSYSTEM  EVALUATION  OPTION  SELECTED) 

GO  TO  8203 

8202  READ(1,*) DSMAX, DVMAX 
WRITE (6,931) DSMAX , DVMAX 

931  FORMAT (/2 0X3 5HSYSTEM  DESIGN  OPTION  SELECTED  WITH: /35X22HMAX  BATTER 
CY  DISCHARGE**,  F5. 1, 2H  (,/3  5X20HMAX  SPEED  DEVIATION-,  F5 . 1 ,  6H  KM/HR) 

8203  DEV-0. 


DIS=0. 

83  W1=0. 

V(1)=0. 

VX(1 )=0 . 

W=0 . 

WB=0 . 

WH=0 . 

D=0 . 

KK=1 

PMAX=PRATM* I PM AX 
MM=5 . 75**ALOG10 (PRATM) *2 1 . 5+HM 
M= ( MM+BM+FM ) / . 74+ . 8  5  *  PM 
AHT=0 . 

IF(DIS. EQ. 0. )G0  TO  8302 

8301  FA=(M/M1 )  ** . 5*FA 

8302  M1=M 
DIS=0. 

DO  160  11=1, NNNDATA 
WRITE (6, 924) 

924  FORMAT (1H1, '  TIME  VELO-SCH  VELO-ACT  DEVIATION  F-ROLL  F-AIR 

C  F-ACCEL  INPUT-POWER  GEAR  SHAFT-RPM  CONTROLLER  MTR  ARM' 
C  '  BAT  DISCHARGE ' ) 

WRITE (6, 9241) 

9241  FORMAT (IX, 1  (SEC)  (KM/HR)  (KM/HR)  (KM/HR)  (NT)  (NT) 

C(NT)  (MTR)  (HYB)  SELECT  (MTR) (HYB)  STEP/FIELD  VOLTS  AMPS 
C  (%)') 

DO  160  111=1,60 

1=111+11*60-59 

IF(VX(I) .GE.99.9JGO  TO  170 

*  85  VBAT=(73.1-4.32*DIS*DIS)*VNOM/72. 

V(I)=VX(I) 

KKK=0 

90  W=V  (I)*. 2778 
DELTAV=. 02 
A=0 . 

VLT=0 . 

MS=0 . 

MSH=0 . 

MMM=0 

*  PRQDH=0 . 

*  PAVL=0 . 

MCF-1 . 02 
WA=(Wl+VV)/2. 

IF( IS . EQ. 0 )G0  TO  92 

91  CALL  MTRSPD(VVA, RR, AR, GR, UPSHFT, DNSHFT, IS, KK, MS) 

MCF-MCF+ . 0002  94  * ( AR*GR ( KK ) /RR ) *  * 2 . 

92  FR= ( CR+ . 000000209*VVA**2 . 8) *M*9 . 81 
FAIR-VVA* *2 . * 1 76 . 4*FA* CD/ ( 2 73 . +TMP ) 

FACC=M*  (W-VV1 )  /T*MCF 
F-FR+FAIR+FACC 

P»F*WA/1000. 

IF (  WA .  EQ .  0 .  ) GO  TO  153 
9002  PSPD=V (I )/VMAX*100. 

PP=P/PMAX*100. 

IF ( PP . EQ . 0 . ) GO  TO  9202 


9201  ETAA=95.8-( . 01217*PSPD*PSPD+. 8879*PSPD+4 . 261 ) /ABS (PP) 

9202  IF( IS . NE. 0)GO  TO  9204 

9203  MS=BSPD 

CALL  TRNSEFF(  P,  ETAA,  PRATT,  WA,  RR,  AR,  BSPD,  TSTP,  ITSTP,  LTAT,  ETA) 

GO  TO  100 

9204  IF( PP. EQ. 0. )GO  TO  153 

9205  PSPDT»PSPD*GR(KK)/GR(IS) 

IF ( GR ( KK ) . GE . 1 . ) GO  TO  94 

93  EMAX*EXP ( . 03* ( 1 . -1 . / GR(KK) ) ) 

GO  TO  95 

94  EMAX=EXP( . 03* ( 1 . -GR (KK) ) ) 

95  ETA=EMAX*(99.-( . 004005*PSPDT*PSPDT+. 1849*PSPDT-1 . 565 ) /ABS (PP ) ) 

100  IF ( PP . LE . 0 . ) GO  TO  102 

101  PRQD«P/ETAA/ ETA* 10000. 

GO  TO  (110, 1101, 1101)CTYP 

102  PAVL=P*ETAA* ETA/ 10000. 

GO  TO  151 

*  110  CALL  VSTEP(PRQD, MS, RARM, KMOT, VBAT, ICSTP, CSTP, A, PAVL, AMAX, PRATM, 

*  CMSTP, IMSTP, ETAM , MMM , RBAT , VLT, IHYB , HR, ETAH, PRQDH , V ( I ) , VH, MTYP, RFLD, 

*  CCOEF1, COEF2, COEF3 ) 

GO  TO  111 

*  1101  CALL  VCHOP(MTYP, RARM, RFLD, KMOT, AMAX, PRATM, IMSTP, MSTP, ETAM, RBAT, VBA 

*  CT , IHYB , HR , ETAH , VH , PRQD , MS , V ( I ) , A, VLT , PAVL , PRQDH , COEF 1 , COEF  2 , COEF  3 ) 

111  IF ( ABS ( ( PAVL-PRQD ) /M/V ( I ) ) .GT. . 0001 )DELTAV=. 5 

112  IF ( KKK . EQ . 2 )GO  TO  121 

118  IF ( PAVL. LE. PRQD) GO  TO  120 

119  V ( I )=V ( I ) +DELTAV 
KKK=1 

GO  TO  90 

120  IF(KKK.LE.l)GO  TO  121 

121  IF ( PAVL. GE. PRQD) GO  TO  152 

122  V(I )=V ( I )— DELTAV 
KKK=*2 

GO  TO  90 

151  WBs*WB+T/3600.  *PAVL 
GO  TO  153 

152  W»4f+T/3600.*PAVL 
WH»WH+T/ 3600.*  PRQDH 

153  IT=I*T-T 
W1=VV 

PRQDE“P AVL -PRQDH 

IF ( V ( I ) . GE . VH . AND . IHYB . EQ. 1 )MSH=MS*HR 
D»D+VVA*T/1000. 

DEV«V(I)-VX(I) 

WRITE(6, 918)IT, VX(I ),V(I), DEV, FR, FAIR, FACC,PRQDE, PRQDH, KK, MS, MSH,M 
CMM, FLD, VLT, A, DIS 

918  FORMAT (1X1 4, F9. 1,2F10.1,F10.0, 2F8. 0, F7. 1, F6. 1, 15, F8. 0, F6. 0, 14, F7. 0 
C,F7.1,F6.0,F12.4) 

IF ( MMM . EQ . 0 ) MMM=1 

155  CALL  BATDIS(TMP, A, AHT, T, IT, DIS, CSTP (MMM) , VNOM, BM) 

IF( IPROG. EQ. 0 )GO  TO  160 
IF (DIS. LE. DSMAX)GO  TO  1553 

1551  BM-BM+100. 

DO  1552  LL-1, ICSTP 

1552  RBAT(LL)*RBAT1(LL)*BM1/BM 


noon 


PRINT*,'  TOTAL  TIME*  ' , IT, '  SECS' 

WRITE (6,928) BM 

928  FORMAT (/20X92H .  CYCLE  TERMINATED  DUE  TO  INSUFFICIENT  BATTERY  C 

CAPACITY  TO  MEET  RANGE  REQUIREMENTS  . ,/35X24HNEXT  ITERATION  WIL 

CL  USE: , / 40X14HBATTERY  MASS*  ,F5.0) 

GO  TO  83 

1553  IF(V(I) .GE.87.  .OR. IT. LE. 120)GO  TO  160 

156  PRATM =PRATM+1 . 

R=PRATM/ (PRATM-1 . ) 

BM=BM*R*  * . 5 

DO  157  L=1 , ICSTP 

157  RBAT(L)*RBAT1(L)*BM1/BM 
RFLD=RFLD/R 
RARM=RARM/R 
KMOT=KMOT/R 
AMAX=AMAX*R 

PRINT*,'  TOTAL  TIME=  *,IT, ’  SECS' 

WRITE (6,925) PRATM , BM 

925  FORMAT ( 2 5X81H ' ' ' ' '  CYCLE  TERMINATED  DUE  TO  INSUFFICIENT  POWER  TO  M 
CEET  POWER  REQUIREMENTS  /3 5X24HNEXT  ITERATION  WILL  USE:/40X13 

CHMOTOR  POWER=  , F4. 1/40X14HBATTERY  MASS=  ,F5.0) 

GO  TO  83 
160  CONTINUE 

170  PRINT*,'  TOTAL  CYCLE  TIME  (MIN)  =  ' 

TIME=T*(IT-1 )/60. 

WRITE (6, 919)TIME 
919  FORMAT (F10. 4 ) 

PRINT*,'  TOTAL  CYCLE  DISTANCE  (KM)  =' 

WRITE ( 6 , 91 9 ) D 
WM=W— WH 

PRINT*,'  TOTAL  MOTOR  INPUT  ENERGY  (KW-HR)  =  ' 

WRITE (6,919) WM 

PRINT*,'  TOTAL  ENGINE  INPUT  ENERGY  (KW-HR)  =' 

WRITE (6, 91 9 )WH 

PRINT*,'  TOTAL  REGENERATION  ENERGY  AVAILABLE  (KW-HR)  =  ' 

WRITE (6,919) WB 

PRINT*,'  BATTERY  DISCHARGE  AT  END  OF  CYCLE  (()  =’ 

WRITE (6, 919)DIS 
GAS=WH* . 56*FCOST/D 
ELECT=4fM*ECOST/  .  6/D 

PRINT*,'  AVG  CYCLE  ELECTRICAL  COST  ($/KM)  =' 

WRITE (6, 91 9) ELECT 

PRINT*,'  AVG  PETROLEUM  FUEL  COST  ($/KM)  =' 

WRITE (6, 919)GAS 

IF ( I PROG . EQ . 0 ) GO  TO  172 

PRINT*,'  FINAL  VEHICLE  PARAMETERS  ARE:' 

PRINT*,'  MOTOR  POWER*  PRATM,'  BATTERY  MASS*  ' , BM 

PRINT*,'  TOTAL  MASS*  ' ,M, '  FRONTAL  AREA*  ' ,FA 

172  END 

SUBROUTINE  MTRSPD(VVA, RR, AR, GR, UPSHFT, DNSHFT, IS, K, RMS) 

DIMENSION  GR (10) 

THIS  SUBROUTINE  COMPUTES  MOTOR  SPEED  FOR  MULTI-SPEED  TRANSMISSIONS 
AND  DRIVING  CYCLE  REQUIREMENTS. 
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100  RMS=VVA*9.549/RR*AR*GR(K) 

IF  ( RMS . LE . UPSHFT )G0  TO  103 

101  IF(K.GE. IS)GO  TO  106 

102  K=K+1 

GO  TO  100 

103  IF  ( RMS . GE . DNSHFT )GO  TO  106 

104  IF(K.LE.l)GO  TO  106 

105  K=K-1 

GO  TO  100 

106  RETURN 
END 

SUBROUTINE  VCHOP(MTYP, RARM, RFLD, RKMOT, AMAX, PRATM, IMSTP, RMSTP, ETAM, 
CRBAT , VB AT , IHYB , HR, ETAH, VH, P, RMS, V, A,VLT, PAVL, PRQDH, COEF1 , COEF2 , COE 
CF3 ) 

DIMENSION  RMSTP ( 10 ) , ETAM( 10, 30 ) 

THIS  SUBROUTINE  CALCULATES  MOTOR  POWER  DEVELOPED  AND  CURRENT 
REQUIRED  FOR  SCR  CHOPPER  VOLTAGE  CONTROLLERS  WITH  DC  MOTORS.  IT 
ASSUMES  1.5  VOLTS  JUNCTION  LOSS  AND  3(  COPPER  LOSS. 

KK=0 

C=RMS*RKMOT 
PAVLH=0 . 

PRQD=P*1000. 

IF( IHYB . NE. 1 )GO  TO  423 
IF ( PRQD . LE . 0 . ) GO  TO  423 
IF(V.LT.VH)GO  TO  423 
*414  HRS=RMS*HR 

*  TORK=«COEFl4COEF2*HRS+COEF3*HRS**2 

*  PA VLH=(TORK*HRS*. 7457/5252. ) *ETAH/100. 

*  IF ( PAVLH . GE . P )GO  TO  421 

*  PRQDH«PAVLH 

*  PRQD=PRQD— PRQDH*1 000 

*  GO  TO  423 
*421  PRQDH*P 

*  PAVL*P 

*  RETURN 

423  CALL  MTREFF( PRQD, PRATM, RMS, RMSTP, IMSTP, ETAM, ETA) 

IF ( MTYP . NE . 0 )GO  TO  402 

401  EAIA=PRQD/ETA*100.*C/(RFLD+RARM*RFLD/(RFLD-C) ) 

IF ( EAIA . LT . PRQD ) E AIA=PRQD 
VLT-(EAIA*RARM*RFLD*RFLD/C/ (RFLD-C) )** .5 

VM AX-AMAX  *  RARM  * RFLD/ ( RFLD-C ) 

A-EAIA*RFLD/VLT/ C 
GO  TO  403 

402  RTOT-RARM+RFLD 
EAIA«PRQD*C/ETA/ ( C+RTOT ) 

IF { EAIA . LT . PRQD ) EAIA=PRQD 
A«(EAIA/C)**.5 

VMAX"AMAX* ( RBAT+1 . 0  3  * ( RTOT+C ) )+1.5 
VLT«A*VMAX/AMAX 

403  IF ( VLT . GT . VM AX ) KK=1 
VLT-VMAX 

A-AMAX 

IF(VLT.GT. VBAT)GO  TO  406 


*onoo 


404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

*  416 
417 

* 

* 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

*  521 

* 

* 

523 

501 

5011 

5111 

5112 

5113 


IF(KK.NE.0)GO  TO  407 

PAVL=P 

RETURN 

VLT=VBAT 

IF ( MTYP . NE . 0 ) GO  TO  409 

EAIA=C* VLT*VLT/RARM/ RFLD* (1. -C/RFLD) 

A=EAIA/C/VLT*RFLD 
GO  TO  410 

EAIA* (VLT-1 . 5 ) / ( RBAT+1 . 03  * ( RTOT+C ) ) **2  .  *  C 
A=(EAIA/C)**.5 

CALL  MTREFF ( E AIA , PRATM, RMS, RMSTP, IMSTP, ETAM,  ETA) 

IF ( MTYP . NE . 0 )GO  TO  412 

PAVL=ETA* (EAIA+A* A*RARM+VLT*VLT/ RFLD) *1 00000 . 

GO  TO  413 

PAVL=ETA* (EAIA+A* A* RTOT ) /I 00000 . 

IF ( PAVL . GT . EAIA/ 1000.) PAVL=EAIA/ 1 000 . 

PAVL=PAVLH+PAVL 

RETURN 

END 

SUBROUTINE  VSTEP(P, RMS, RARM, RKMOT, VBAT, ICSTP, CSTP, A, PAVL, AMAX, PRAT 
CM, RMSTP, IMSTP, ETAM, M, RBAT, VLT, IHYB , HR, ETAH, PRQDH, V, VH, MTYP, RFLD, 
CCOEF1 , COEF2, COEF3 ) 

DIMENSION  CSTP (5 ) , RMSTP ( 10 ) , ETAM( 10, 30) ,RBAT(5 ) 

THIS  SUBROUTINE  CALCULATES  MOTOR  POWER  DEVELOPED  AND  CURRENT 
REQUIRED  FOR  FINITE  STEP  VOLTAGE  CONTROLLERS  WITH  DC  MOTORS. 

PRQD=P 
PAVLH=0 . 

IF( IHYB .NE. 1 .OR. V. LT. VH)GO  TO  523 
IF ( PRQD . LE . 0 . ) GO  TO  523 
HRS=RMS*HR 

TORK=COEFl+COEF2*HRS+COEF3*HRS**2 
PAVLH= (TORK*HRS*  .  7457/5252.  )  *ETAH/100. 

IF ( PAVLH . GE . P )GO  TO  521 
PRQD=P-PAVLH 
PRQDH=PAVLH 
GO  TO  523 
PRQDH-P 
PAVL=P 
RETURN 
MM-0 

IF(MTYP. EQ.0. AND. RMS. LT. 1000. )RMS=1000. 

M»M+1 

C=RMS*RKMOT 

DV»A*RBAT(M) 

VLT-CSTP(M) *VBAT-DV 
IF(MTYP. NE . 0)GO  TO  5112 
EAIA»C*VLT*VLT/ RARM/ RFLD* (1 . -C/RFLD) 

A»EAIA*RFLD/VLT/ C 
GO  TO  5113 
RTOT-RFLD+RARM 

EAIA-( (VLT-VLT*C/( RTOT+C) ) /RTOT) **2 . *C 

A-(EAIA/C)**.5 

DV2-A*RBAT(M) 
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IF( ABS (DV2-DV) .LE. 0. 3 )GO  TO  5013 

5012  DV=DV2 

GO  TO  5011 

5013  CALL  MTREFF(EAIA, PRATM, RMS, RMSTP, IMSTP, ETAM, ETA) 

IF( MTYP. NE . 0 )GO  TO  5212 

5211  PAVL -ETA* (EAIA+A*A*RARM+VLT*VLT/RFLD) /100000. 

GO  TO  5213 

5212  PAVL=ETA* (EAIA+A*  *2 . *RTOT ) /l 00000 . 

5213  IF( PAVL. GT . EAIA/1000 . ) PAVL=EAIA/1000 . 

5015  IF( MM. NE. 0)GO  TO  508 

502  IF( M.GT . 1 . AND. A. GT . AMAX )GO  TO  506 

504  IF(PAVL.GE. (PRQD* .8) )GO  TO  5061 

505  IF ( M . LT . ICSTP )GO  TO  501 
5061  IF ( PAVL. GE. PRQD) GO  TO  507 

IF(M.NE. ICSTP )GO  TO  508 

506  M=M-2 
MM=1 

GO  TO  501 

*  508  PAVL =PAVL  +P A VLH 
RETURN 

507  PAVL=P 
RETURN 
END 

SUBROUTINE  MTREFF( EAIA, PRATM, RMS, RMSTP, IMSTP, ETAM, ETA) 

DIMENSION  ETAM( 10, 30) , RMSTP(10) 

THIS  SUBROUTINE  PERFORMS  A  BILINEAR  INTERPOLATION  FOR  MOTOR 
EFFICIENCY(ETAM)  AS  INPUTTED  FOR  DISCRETE  VALUES  OF  MOTOR  SPEED 
(RMSTP)  AND  PERCENT  OF  RATED  POWER ( PP ) . 

PP-EAIA/ PRATM/ 1 0 . 

K=1 
KK=1 

600  IF(RMS.LE. RMSTP (K ) )GO  TO  601 

601  IF (K.GE. IMSTP )GO  TO  603 

602  K=*K+1 

GO  TO  600 

603  IF(PP.LE. (KK*10. ) )GO  TO  605 
6031  IF(KK.GE. 30)GO  TO  605 

604  KK=KK+1 
GO  TO  603 

605  IF(K. NE . 1 )GO  TO  607 

606  IF(KK.EQ. l)GO  TO  608 
GO  TO  609 

607  CC* ( RMS-RMSTP (K-l ) ) /(RMSTP(K)-RMSTP(K-l ) ) 

ETA2* ( ETAM ( K , KK ) -ETAM ( K -1 , KK ) ) *CC+ETAM ( K-l , KK ) 

IF(KK. EQ. 1 )GO  TO  610 

GO  TO  611 

608  ETA»ETAM(1, 1) 

RETURN 

609  ETA1«ETAM(K,KK-1) 

ETA2«ETAM(K,KK) 

GO  TO  612 

610  ETA-ETA2 
RETURN 


611  ETA1=(ETAM(K, KK-1 )-ETAM(K-l , KK-1 ) ) *CC+ETAM(K-1 , KK-1 ) 

612  ETA= (ETA2— ETA1 ) * (PP— (KK-1 ) *10. ) /10 . +ETA1 
RETURN 

END 

SUBROUTINE  TRNSEFF(  P ,  ETAA,  PRATT,  WA,  RR,  AR,  BSPD,  TSTP,  ITSTP,  ETAT,  ETA 
C) 

DIMENSION  TSTP( 10 ) , ETAT (10, 20 ) 

THIS  SUBROUTINE  PERFORMS  A  BILINEAR  INTERPOLATION  FOR  CVT  EFFICI- 
ENCY(ETAT)  AS  INPUTTED  FOR  DISCRETE  VALUES  OF  SPEED  RATIOS(TSTP) 
AND  PERCENT  RATED  POWER (PP). 

PP=P*ETAA/ PRATT 

DSS=9. 549*VVA/RR*AR 

SR=DSS/BSPD 

K=1 

KK=1 

700  IF(SR.LE.TSTP(K) )GO  TO  703 

701  IF (K.LT. ITSTP )GO  TO  703 

702  K=K+1 

GO  TO  700 

703  IF( PP . LE. (KK*5 . ) )GO  TO  705 
7031  IF ( KK . NE . 20 ) GO  TO  705 

704  KK=KK+1 
GO  TO  703 

705  IF(K.NE. 1 )GO  TO  707 

706  IF ( KK . EQ . 1 ) GO  TO  708 
GO  TO  709 

707  CC=(SR-TSTP(K-1 ) )/(TSTP(K)-TSTP(K-l ) ) 

ETA2* (ETAT (K , KK ) -ETAT (K-l , KK ) ) *CC+ETAT (K-l , KK ) 

IF(KK. EQ. 1 )GO  TO  710 
GO  TO  711 

708  ETA=ETAT (1,1) 

RETURN 

709  ETA1=ETAT(K, KK-1 ) 

ETA2=ETAT(K,KK) 

GO  TO  712 

710  ETA=ETA2 
RETURN 

71 1  ETA1= ( ETAT (K , KK-1 ) -ETAT ( K-l , KK-1 ) ) *CC+ETAT ( K-l , KK-1 ) 

712  ETA*(ETA2-ETA1 ) * ( PP- (KK-1 ) *5 . ) /5 . +ETA1 
RETURN 

END 

SUBROUTINE  BATDIS (TMP, A, AHT, T, IT, DIS, CSTP, VNOM, BM) 

THIS  SUBROUTINE  COMPUTES  TOTAL  BATTERY  DISCHARGE )(  (DIS)  FOR  LEAD 
ACID  BAT  ERIES  OF  GOLF-CART  TYPE  CONSTRUCTION.  IT  IS  TEMPERATURE 
CORRECTED  AND  USES  A  MODIFIED  FRACTIONAL-UTILIZATION  METHOD  WITH 
CORRECTIONS  FOR  CHANGING  CURRENTS  AND  RECUPERATION  PERIODS. 

A=A*CSTP*350./BM*VNOM/72. 

AHT»AHT+A*T/ 3 600 . 

AH«(195 . -. 633*A+. 00091 3*A* A) * ( . 014*TMP+. 65 ) 

AAVG«AHT/IT*3600. 

AHA»(195.-.633*AAVG+. 00091 3*AAVG*AAVG) * ( . 014*TMP+. 65 ) 


IF( A.NE. 0. )G0  TO  802 

801  RA=0 . 

GO  TO  805 

802  IF( A. LE. AAVG )GO  TO  804 

803  RA=AAVG/A 
GO  TO  805 

804  RA=A/AAVG 

805  DIS=(AHT/AH*RA+AHT/AHA* (1 . -RA) ) *100. 
RETURN 

END 
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